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ABSTRACT 


The usefulness of a barotropic prediction model for providing real- 
istic displacements at levels other than 500 mb is suggested to circumvent 
difficulties with current baroclinic models attempting to predict develop- 
ment as well. 

A prediction equation is obtained by making use of concepts from 
Charney's equivalent barotropic model, ‘Arnason's stratified model, and 
Cressman's semi-empirical barotropic model. Further, a stream—-function 
prediction scheme, making use of geostrophic and non-divergent wind rela- 
tionships, is used to make a series of 24- and 48-hour forecasts of 300-mb 
heights. The vorticity equation used is 
~S Kae) + + JOY, Kao +f) - Kae TQ, y)= O 
where K;, Ko, and K3 are empirical constants. Results are presented for 
the effects of adjusting Ky and K3. 

The author wishes to express his appreciation to Professor George J. 
Haltiner of the U. S. Naval Postgraduate School for his guidance and en- 
couragement in this investigation. 

Appreciation is also expressed to the personnel of the U. S. Fleet 
Numerical Weather Facility for their assistance in this project. Special 
thanks are extended to Lieutenants Mildred J. Frawley and Harry &. 
Nichoison and Mr. Milton H. Reese for assistance in programming and many 


helpful suggestions. 
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1. <daimaduction. 

Up to the present time, most attempts to predict the pressure height 
at levels other than 500 mb have been made with baroclinic models. The 
results have been disappointing because of inadequately controlled devel- 
opment mechanisms. Observing the persistent skill of barotropic forecasts 
at 500 mb, 'Arnason [1] suggested the possiblity of devising a barotropic 
prediction model to be applied to other levels, "capable of realistic dis- 
placements of pressure systems...but containing no mechanism for develop- 
ment," 

This investigation introduces an empirically adjusted divergent baro= 
tropic model for application to the 300=-mb surface. Concepts from Charney's 
[2] equivalent barotropic model and 'Arnason's stratified model [1] are 
incorporated in addition to empirical results by Cressman [3] on control 


of very long atmospheric waves. 





2. Background, 

Looking first to the equivalent barotropic model, it is recalled that 
partial baroclinicity is permitted in that wind speed, but not direction, 
is allowed to vary with height at any point. If the integrated mean of 
this wind with respect to pressure is determined, the wind, Y, and vor- 
ticity,  , at any level can be expressed as a scalar factor, A, of this 
mean. For further simplicity, A is considered to be a function of pressure 
alone. Thus 

VENER. ‘ge 2 
Y =ACPyY , C=Ape , (2.1) 
where the cap indicates the vertically integrated mean. Climatologically, 
y is observed to occur at about 600 mb. 
Substituting the expressions (2.1) in the approximate vorticity 


equation 


4 vy (Er) =f. 


and integrating with respect to pressure gives 


ce 
Ss> 


+ UY -9(ME +f) =0. ED 


Qy 
c+ 


Here f is the coriolis parameter, and & qemeena empirically to have a 
value of about 1.25. Multiplying (2.2) by i and substituting rn e = Cand 
RY = ‘V* forms the basis for barotropic forecasts at the equivalent baro=- 
tropic level, (p*), occurring in nature at about 500 mb. 


Nn A 
If, however, the expressions (2.1) are substituted for VY and Z into 


(2.2), the following is obtained: 


a ee K.0+4) =O (2.5) 





/ 
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where Ky = 4. For the wind field (2.1), equation (2.3) applies to any 
level with the Kp appropriate to that level. 
In applying the theory of the stratified model to forecasting, 


'Arnason arrives at an expression for divergence 


9-V =-3 (28 4 9-94) : (eh) 
where b is geopotential, VY is the basic flow, and H is the thickness of 
the stratified layer. The parameter, g', is related to gravity, although 
much smaller, and depends upon static stability and the vertical extent 
of the stratified layer. This expression for divergence is introduced into 


the vorticity equation 


— 4. YW 3 \V "a + "Vv -Y O 75) 
to give 


1 Vv: ~— 6 
Ne 4 Wma — sh ( Se + V-VS) =O, (2.6) 


where }) =f +# 4s absolute vorticity; YJa is defined for our purposes as 
>. | soe 
4 = Kf A { of (2esee 
The prediction scheme consists of introducing a stream function, WY , 


to represent the wind and vorticity 


VekxVy, S2Vy (267) 


where le is the unit vector in the vertical. The stream function is ini- 


tially obtained from the balance equation 


FTA y rd (Wer Way — Vey) + VY VE - 7°G =O - (2.8) 


Next, in order to simplify the procedure, the quasi-geostrophic relation- 





ship 


) rs {J (2.9) 


is introduced into the divergence term. The result of substituting (2.7) 


and (2.9) into (2.6) is 
3 vy + TOV) - BLS + T,W)] =o 
(ve =e + JOY, Na) es: JO, y) =O J (2.10) 


where Y) is the stream function corresponding to the basic flow. Equation 


(2.10) is a prognostic equation in a single variable, YW. 





3. Procedure. 

Application of ‘Arnason's stratified model involves an interpretation 
of the parameters, YJ, g', and H. Also, the introduction of a factor, Kj, 
into the tendency portion of the divergence term, and a factor, K3, into 
the advection portion of the divergence term permits adjustment of the mag- 
nituce of these terms for various levels. 

The basic flow for this prediction model is taken as a spaced-averaged 
flow at 300 mb, employing a 15=pass, 4-point smoother over the finite-dif- 
ference grid used for computations (grid distance equal to 381 km at 60N 


latitude): 


Wo = ma + Wa + Ws + Wu ) (3.1) 


The subscript zero denotes the central point of a diamond array and the 
subscripts one through four denote the surrounding four grid points, 

In arriving at an interpretation of g', reference is made to Cressman's 
[3] treatment of a model including a dependence upon static stability. He 
introduced a stability parameter, p, which as applied here amounts to a 
substitution of 

_— Of 

In a series of 500-mb forecasts Cressman used various values of , in his 
divergent barotropic model. It was concluded that = affords the de= 
Sired control on the ultra-long waves, and gives better forecasts than a 
former Joint Numerical Weather Prediction Unit forecast model. Use has 
been made of this value by the substitution of (3.2) into (2.10). 

In addition, the mean value 9 25 has been used for 3p (similar to 


Fen 
Cressman), taking %,,,.as the U. S. standard atmosphere 500-mb height, 





and ae f at 45N latitude. 


With these last adjustments, (2.10) becomes 


(78 — keel) SY + TOG S+#) -— KARI TCY) =O. 9) 


In finite differerce form, (3.3) becomes 


(VA KebI ae + £0, bf) — Kyfeet NP V)=0 , (3.4) 


where 7%, A , and are the finite-difference equivalents of the Lapla- 
cian, differential, and Jacobian operators, respectively; and gd is the 
sea-level distance between grid points. Forecasts are made in one-hour 
time steps, using a forward difference for the first step and centered 
differences for successive steps. 

Equation (3.4) was programmed for the Control Data Corporation Model 
1604 computer using U. S. Fleet Numerical Weather Facility's 1977-point 
octagonal grid, which covers the Northern Hemisphere from pole to about 
LON. 

Step one in the forecasting procedure consists of solving the balance 
equation for an initial stream function, using analyzed D=values (deviation 
in height from standard). A technique formulated and programmed for the 
CDC 1604 by 'Arnason and Reese at the Fleet Numerical Weather Facility was 
used. It was found that for most of the 300-mb data available for this 
investigation, this method would not converge to the cut-off criterion 
(all residuals less than 0.5 feet in the relaxation procedure as normally 
used at lower levels). In order to proceed with the work, an arbitrary 
cut-off was introduced if non-converrence became apparent. This intro- 
duced errors at a few points of the grid due to residuals which were now 


of the order of three to five feet. These were considered to be insignif- 





icant with respect to the prediction scheme, 

Step two was solution of equation (3.4), and this was repeated on 
each day's data using various values of the factors Ko and K3- The factor 
K, was taken equal to one and not adjusted in this investigation. Values 
of Ky and K3 were chosen as work progressed so as to reduce forecast er= 
rors to a minimum, 

Step three consisted of verifying the forecasts against the prepared 
analyses by a simple point-by-point difference, and calculation of the 


pillow and root-mean-square error (RMS#): 

5 ks a a 94 

Pillov = 4 (A-B)n > RMSE =[4>(A-B-Pillow), | (3.5) 
nz] nei 


x = 1977 points, A = forecast map, B = analysis map. 





4. Results. 

Table 4.1 lists the results of adjusting the factor Kz with Kp held 
equal to one. Numbers to the left of each column are pillow while the 
numbers to the right are RMSE. Upper numbers of each paired set are for 
24-hour forecasts and lower numbers are for 48-hour forecasts. 

Runs were not made for all included values of K3 but minimum errors 
were sought. The numbers in parentheses are values taken from faired 
curves through the values that were run. Fig. 4.1 shows the curves for 
24-hour forecasts; those for 48-hour forecasts would be similar. The 
trends of errors for the successive values of K3 were observed to be quite 
regular, thus justifying an approximation of intermediate values where 
necessary for the averages over all days. In cases where actual reversal 
of error trends was not observed, the fact that differences between suc= 
cessive values of K3 were reduced to zero or one foot was taken as an in- 
dication that reversal was imminent, and adjustments of K3 were terminated. 

In one case (17 January) a forecast was made using a positive value 
of I'3, showing an error much larger than the optimum for this day. Other 
work on this term not presented here further confirms the large errors for 
values of K3, including K3 equal zero (exclusion of the term), beyond the 
range indicated herein. 

The column entitled "Balancing errors" indicates the differences be- 
tween (a) the height fields obtained after balancing to the stream func- 
tion and inverting again to height, and (b) the initial height field. 

Table 4.2 lists the results of adjusting the value of Kp. Since K3 
was first adjusted, optimum values of this factor were used in adjusting 
Ko. 


It will be noted from table 4.1 that a considerable portion of the 


8 





forecast errors is apparently due to the stream-function conversion, On 
13 January, the conversion error amounts to 41% of the total forecast er- 
ror, while on other days the ratios of errors range from 17% to 25%. 

The trends of errors for varying K3 all indicate a minimum error in 
the range -0.5 to -1.75 and also indicate the advantage of including 
'Arnason's divergence term, From extrapolation of the curves in fig. 4.-, 
errors for K3 equal zero are estimated to average about 20 feet above the 
least average for 24-hour forecasts. 

The use of the factor Ky does not seem to justify any deviation from 
unity, with improvements ranging from zero to only three feet for 24-hour 


forecasts, 
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Figs. 4.2 through 4.4 depict the 13 January analyzed height field, 
the 14 January analyzed height field, and the 13 January, 24-hour forecast 
height field verifying on 14 January, respectively. In fig. 4.4 the ini- 
tial (13 January) and observed (14 January) characteristics of certain of 
the systems are superimposed upon the forecast (for 14 January) field to 
indicate contrasting features. Initial features are dotted, and observed 
are dashed. 

The following characteristics may be noted: 

(1) In the western Pacific, the forecast of the short-wave trough 
was very much in agreement with its actual movement and intensity. No 
cyclogenesis occurred in this system and forecast errors were near zero 
in the area of the trough. 

(2) The eastern Pacific low occurring on 13 January almost completely 
diss .pated in 24 hcurs, The forecast position of this low is nearly that 
of the 14 January trough identifying the system. Height errors in excess 
of 700 feet are due to cyclolysis, which of course is not provided for in 
the model. 

(3) In the Alaskan ridge, cyclogenesis occurred from 13 to 14 Janu- 
ary in the NW portion, accounting for errors up to 600 feet. The ridge, 
however, was well forecast. 

(4) The trough off the Atlantic coast of the U. S. was forecast to 
move in agreement with observation, but cyclogenesis again resulted in er- 
rors of about 400 feet. 

(5) The low initially over the central Mediterranean moved into a 
cut-off position south of a blocking high. The forecast movement was in 
fair agreement, although a bit short, but accompanying cyclogenesis raised 


errors to about 600 feet. The forecast position of the blocking high was 


abs 





good, but again, anticyclogenesis in the high and the ridge to the north- 
east introduced forecast errors of 400 feet and 600 feet, respectively, 
“isewhere on the map, the significant errors seem generally to be due 
to development processes and not to movements of the systems. These ob- 
servations were corroborated by investigation of another sequence of maps 


produced in the investigation, 
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5. Conclusions. 

It is noted from the averaged errors in table 4.1 that there are only 
limited differences in results over the range of Ky values indicated. It 
was apparent, however, that sizable differences occur as K3 is extended 
beyond this range, and of particular interest is the advantage over exclu- 
sion of the term (K3 equals zero). These results indicate an optimum val- 
ue for K3 of about -1.0. 

Even though varying Kp did not yield significant differences in this 
investigation, it is possible that a greater number of cases might do so, 
particularly if some other season is considered. Only data from December, 
1961 and January, 1962 were used, 

Taking note of the large errors encountered in converting to a stream 
function, it is possible that a geostrophic model, not requiring such con- 
versions, might reduce total errors somewhat. It is hoped that this in- 
vestigation can be accomplished in the near future. 

The model developed in this paper seems to fulfill its intended pur- 
pose, that is, forecasting realistic displacements unaffected by uncon- 
trolled developments. with an understanding of its limitations, a fore- 
caster should have here a useful tool, when supplemented by techniques for 
development prediction, to forecast pressure systems at the 300-mb level 


more easily and accurately. 
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